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Abstract 

We report preliminary results on the analysis of the three-body T(10860) — > BBtt, Y(10860) — > 
[BB* + c.c.]7T and T(10860) B*B*tt decays including an observation of the T(10860) 
^(10610)7^ -». [BB* + c.c.] ± ^ and T(10860) -»• ^(10650)^ -». [B+B^if? decays as inter- 
mediate channels. We measure branching fractions of the three-body decays to be ,6(Y(10860) — > 
[BB* + c.c.^vrF) = (28.3 ± 2.9 ± 4.6) x 10~ 3 and £(1(10860) -»• [5* J B*] ± vr T ) = (14.1 ± 1.9 ± 
2.4) x 10~ 3 and set 90% C.L. upper limit £(T(10860) -»• [BB]^) < 4.0 x 10~ 3 . We also report 
results on the amplitude analysis of the three-body Y(10860) — > Y(nS)ir + Tr~ , n = 1,2,3 decays 
and the analysis of the internal structure of the three-body T(10860) — > hb(mP)ir + ir~ , m = 1,2 
decays. The results are based on a 121.4 fb _1 data sample collected with the Belle detector at a 
center-of-mass energy near the T( 10860). 

PACS numbers: 14.40.Pq, 13.25. Gv, 12.39.Pu 
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INTRODUCTION 



Two new charged bottomonium-like resonances, Zf,(10610) and Zb(10650), have recently 
been observed by the Belle Collaboration in decays of T( 10860) to five different final states: 
T(nS)7T + 7T~ , n = 1,2,3 and hb(mP)7r + ir~ , m = 1,2 [3, 0]. The analysis of the quark 
composition of the initial and final states allows to assert that these hadronic objects are 
the first examples of states of an exotic nature: Zj, should be comprised of (at least) four 
quarks. Several models have been proposed to describe the internal structure of these states. 
One suggests [H that Z&(10610) and Zj,(10650) states might be a loosely bound BB* and 
B*B* systems, respectively. The proximity of the Zb(10610) and Z&(10650) masses to those 
of the sum of the B and B* mesons and the sum of the two B* mesons, respectively, supports 
this hypothesis. In this case, it would be natural to expect that the Zfe(10610) and Z&(10650) 
states decay respectively to BB* and B*B* final states with substantial rates. 

Evidence for the three-body T(10860) — > BB*n decay has been previously reported by 
Belle in Ref. |4| with a data sample of 23.6 fb _1 . In this analysis we use 121.4 fb _1 of data 
accumulated by the Belle detector at a center-of-mass (cm.) energy near the T(10860) to 
studythree-body T(10860) -»■ [B& B^n* decays and to search for T(10860) -»■ Zfir* -> 
[B^B*]^ decays. 

Note that we reconstruct only three-body B^B^rr combinations with a charged primary 
pion. For brevity, we adopt the following notations: the sum of B + B°n~ and B~B°7r + 
final states is referred to as BBn; the combination of B + B*°7T~, B~B*°tt + , B B*~tt + and 
B°B* + n~ final states is referred to as BB*n and the sum of B* + B*°tt~ and B*~B*°tt + final 
states is denoted as B*B*n. 

THE BELLE DETECTOR 

The Belle detector |5j is located at the single interaction point of KEKB |6J, an asymmetric 
energy double storage ring collider. The detector is a large-solid-angle magnetic spectrometer 
based on a 1.5 T superconducting solenoid magnet. Charged particle tracking is provided 
by a four-layer silicon vertex detector and a 50-layer central drift chamber (CDC) that 
surround the interaction point. The charged particle acceptance covers laboratory polar 
angles between 9 = 17° and 150°, corresponding to about 92% of the total solid angle in the 
cm. frame. 

Charged hadron identification is provided by dE/dx measurements in the CDC, an array 
of 1188 aerogel Cherenkov counters (ACC), and a barrel-like array of 128 time-of-flight 
scintillation counters (TOF); information from the three subdetectors is combined to form 
a single likelihood ratio, which is then used in kaon and pion selection. Electromagnetic 
showering particles are detected in an array of 8736 CsI(Tl) crystals (ECL) that covers the 
same solid angle as the charged particle tracking system. Electron identification in Belle 
is based on a combination of dE/dx measurements in the CDC, the response of the ACC, 
and the position, shape and total energy deposition (i.e., E/p) of the shower detected in 
the ECL. The magnetic field is returned via an iron yoke that is instrumented to detect 
muons and K\ mesons. We use a GEANT-based Monte Carlo (MC) simulation to model 
the response of the detector and determine its acceptance [7(. 
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BACKGROUND SUPPRESSION 



The dominant background comes from e + e _ — > cc continuum events where real D mesons 
produced in e + e~ annihilation are combined with random particles to form a B candidate. 
This type of background is suppressed using variables that characterize the event topology. 
Since the momenta of two B mesons produced from a three-body T(10860) decay are low in 
the cm. frame, their decay products are essentially uncorrelated and the event tends to be 
spherical. In contrast, hadrons from continuum events tend to exhibit a two-jet structure. 
We use ^thr, the angle between the thrust axis of the B candidate and that of the rest of 
the event, to discriminate between the two cases. The distribution of | cos# t hr| is strongly 
peaked near | cos# t hr| = 1-0 for cc events and is nearly flat for B^B^ir events. We require 
|cos#thr| < 0.80 for B — > D^n final states; this eliminates about 81% of the continuum 
background and retains 73% of the signal events. 

Another significant background comes from events with radiative return to a lower 
mass T(4S') state with a subsequent T(4S) — > BB decay. Momenta of B mesons pro- 
duced in this process fall in the same region as those for B mesons from the three-body 
T(10860) B^B^n decays. B mesons originating from the two-body T(10860) B*B*, 
BB*, BB decays produce peaks around P(B) = 1.07 GeV/c, 1.18 GeV/c and 1.28 GeV/c, 
respectively. Momenta of B mesons from three-body T(10860) B^B^ix decays are less 
than 0.9 GeV/c. 

EVENT RECONSTRUCTION 

Charged tracks are selected with a set of track quality requirements based on the number 
of CDC hits and on the distances of closest approach to the interaction point (IP). Tracks 
originated from B candidate are required to have momenta transverse to the beam be greater 
than 0.1 GeV/c to reduce the low momentum combinatorial background. For charged kaon 
identification, we impose a requirement on the particle identification variable, which has 86% 
efficiency and a 7% fake rate from misidentified pions. Charged tracks that are positively 
identified as electrons or protons are excluded. Since the muon identification efficiency 
and fake rate vary significantly with the track momentum, we do not veto muons to avoid 
additional systematic errors. 

Photons from neutral pions are required to produce clusters in the ECL with an energy 
deposition of greater than 50 MeV and not be associated with charged tracks. The invariant 
mass of the two-photon combination is required to be within 12 MeV/c 2 of the nominal 7r° 
mass. The K*° is reconstructed in the K*° — > K + ir~ mode, the invariant mass of the K*° 
candidate is required to be within 70 MeV/c 2 of the nominal K*° mass. The invariant mass 
of the J/ip — > candidates is required to satisfy \M(^i + ^r ) — Mj/^\ < 30 MeV/c 2 , where 

Mj/^p is the nominal mass of the J/ip meson. Neutral (charged) D mesons originating from 
B decays are reconstructed in the D° — > K + n~ and D° — > K + 7t + 7t~tt~ (D~ — > K + 7t~tt~) 
modes. Those originating from D*~ decays are also reconstructed in the D° — > K + tt~7t° 
mode. To identify D*~ candidates we require \M{DH~) - M(D°) - AM D \ < 2 MeV/c 2 , 
where M(D°) and M(D°ir~) are the reconstructed masses of the D° candidate and D°ir~ 
system, respectively, and AMo = Md* — Mb- 

B decays are reconstructed in the following channels: B + — > J/ipK + , B + — > D°tt + , 
B° ->■ J/ijjK* , B° ->■ D~7r + , B° -)• D*~n + . We identify B candidates by their invariant 
mass M(B) and momentum P{B). M(B) and P{B) distributions for B candidates in 
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FIG. 1: (a) Invariant mass, (b) momentum and (c) M r (Bir) distributions for selected B candidates 
in data. Hatched histograms in (b) and (c) show distributions for events in M(B) sidebands. 

data are shown in Figs. [H(a) andQ](b). We require M(B) to be within 30 to 40 MeV/c 2 
(depending on the B decay mode) of the nominal B mass. Mass sidebands are defined as 
50 MeV/c 2 < \M(B) - M B \ < 80 MeV/c 2 . 

Reconstructed B + or B° candidates are then combined with a tt~ candidate and a recoil 
mass to the Bit combination, M r (Bir) } is calculated as M r (Bir) = a/ E^ ms — Pg n , where 
E cms is the cm. energy and Pb-k is the measured three-momentum of the Bit combination. 
Signal T(10860) — > BB*tt events produce a narrow peak in the M r (B7i) spectrum around 
the nominal B* mass, while T(10860) — > B*B*ir events produce a peak shifted to higher 
mass by about 45 MeV/c 2 due to a missed photon from the B* — > Bj decay. It is important 
to note here that, according to signal MC, BB*tt events where the reconstructed B is the 
one from B* produce a peak in the M r (B7r) distribution at virtually the same position 
as BB*tc events, where the reconstructed B is the prompt one. To remove a correlation 
between M r (Bir) and M(B) and to improve the resolution, we use M r (BTc) + M(B) — Mb 
instead of M r (B7c). The M r (Bir) + M(B) — Mb distribution for experimental data is shown 
in Figure dKc), where clear peaks are visible in the BB*ir and B*B*7r signal regions. 

To determine the distribution of background events we combine a reconstructed B candi- 
date with pions of the wrong charge. The M r (B7c) + M(B) — Mb distribution for wrong-sign 
combinations is shown as a hatched histogram in Fig. [T](c). While wrong-sign Bit combi- 
nations reproduce the shape of the combinatorial background very well, the amount of the 
background is underestimated by about 18%. To correct for this effect, we introduce a scale 
factor of 1.18 for wrong-sign combinations; the M r (B7r) distributions shown in Fig. [T](c) 
include this correction factor. The resolution of the signal peaks in Fig. [T](c) is fixed at 
6.1 MeV/c 2 as determined from signal MC. 

ANALYSIS OF T(10860) -> [B^B*]*^ 

The fit to the M r (B7i) + M(B) — Mb distribution for signal events shown in Fig. [T](c) 
yields N B b-k = 1±14, N B b*-k = 184 ±19 and N b *b*-k = 82 ±11 signal events. The statistical 
significance of the observed BB*ir and B*B*ir signal is 9.3a and 5.7a, respectively. The 
statistical significance here is calculated as \f—2 ln(£ /£sig), where £ sig and C denote the 
likelihood values obtained with the nominal fit and with the signal yield fixed at zero, 
respectively. 

For the subsequent analysis of the internal structures of the three-body decays, we require 
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FIG. 2: PDFs for reconstruction efficiency for the (a) T(10860) -)• BB*ir and (b) T(10860) ->■ 
B*B*ir signals. Points with error bars show signal MC with the phase-space distribution of signal 
events, the solid line is the result of the fit, and the dashed line is the phase space with uniform 
efficiency. 

\{M r (B) + M(B) - M B ) - M B *\ < 0.015 GeV/c 2 to select T(10860) ->■ BB*it events and 
\(M r (B) + M(B) - M B ) - (Mb* + £ 7 )| < 0.015 GeV/c 2 , where £ 7 = 0.049 GeV [8], to 
select T(10860) — > B*B*ix events. For selected B^B^n candidate events, we calculate the 
mass recoiling against the charged pion: M r (n) = ^/ E^ ms — P 2 , where Pb-k is the measured 
three-momentum of the charged pion. The M r (ir) distributions for signal T(10860) — > BB*ir 
and T( 10860) — > B*B*n MC events generated with the uniform phase space distribution 
are shown in Fig. [2J To parameterize the M r (7r) dependence of the reconstruction efficiency 
EBB*n( m ), we use the following empirical function: 

E B MBMn( m ) = fl o(l + ai&m + 0-2^n + 03<$m + «4^) X PHSP B (*) bW^w) , (1) 

where m = M r (ii), 5 m = m — m , m = M B (*) + M B * and aj are fit parameters, and 
PHSP(m) is the phase space function for the T(10860) ->• B^B^n decay. To account 
for the instrumental resolution, we smear the efficiency given by Eq. ([T|) with a Gaussian 
function. The resolution a m of the Gaussian is dominated by the cm. energy spread and 
fixed to be a m = 6 MeV/c 2 . The results of the fit are shown in Fig. [2j 

The M r (n) distributions for wrong-sign Bit combinations for events in the BB*n and 
B*B*7i signal regions are shown in Fig. We use the following empirical function to 
parameterize the distribution of background events 

B BWB . v (m) = b e- a5 - x E BMB ^(m), (2) 

where bo and a are fit parameters. As in the case of the fit to the efficiency PDF, the 
background PDF is convolved with a resolution function. Results of fits to sideband events 
are shown in Fig. [31 

The M r (ir) distributions for right-sign Bit combinations in the BB*ir and B*B*n signal 
regions are shown in Fig. HI Excesses of signal events over the expected background levels 
at lower mass edges of the M r (n) spectra are clearly visible for both final states. The 
distribution of signal T( 10860) — > BB*tc events is parameterized with the following model 

S B B*A m ) = (^z b (io6io) + A NR ) x E BB ^(m) } (3) 
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FIG. 3: M r (ir) distribution for wrong-sign Btt combinations for the (a) BB*ir and (b) B*B*ir 
candidate events. Points with error bars are data, the solid line is the result of the fit with a 
function of Eq.([2]). 




FIG. 4: M r (n) distribution for right-sign Bit combinations for (a) Y(10860) — > BB*ir and (b) 
T(10860) — > B*B*ir candidate events. Points with error bars are data, the solid line is the result of 
the fit with the nominal model (see text), the dashed line - fit to pure non-resonant amplitude, the 
dotted line - fit to a single Z^ state plus a non-resonant amplitude, and the dash-dotted - two Z}, 
states and a non-resonant amplitude. The hatched histogram represents background component 
normalized to the estimated number of background events. 

where A^r is the non-resonant amplitude parameterized as a complex constant and the 
Zb(10610) amplitude is a Breit-Wigner function. As a variation of this nominal model, we 
also add a second Breit-Wigner amplitude to account for possible Zfe(10650) — > BB*n decay. 
We also fit the data with only the Zb(10610) channel included in the decay amplitude. The 
results of these fits are shown in Fig. H|a). Two models give about equally good description 
of the data: nominal model and a model with additional non-resonant amplitude. However, 
we select the former one as our nominal model since adding a non-resonant amplitude does 
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TABLE I: Summary of fit results to the M r (7r) distribution for three-body T( 10860) — > BB*ir and 
T(10860) -> B*B*ir decays. 



Mode 


Parameter 


Nominal model 


Model- 1 


Model-2 


Model-3 


Model-4 


BB*ir 




1.08 ±0.12 




0.86 ±0.15 


1.0 


0.73 ±0.17 






0.25 ±0.10 








0.087 ±0.061 




0Z 6 (1O65O) 


-0.93 ±0.34 


— 


— 


— 


0.32 ± 0.23 




JnR 




1.0 


1.37 ±0.28 




1.17 ±0.27 




<Pnr 






0.18 ±0.21 




2.75 ± 1.03 




-log£ 


142 


226 


129 


162 


126 


B*B*7T 


/z b (10650) 


1.0 




0.83 ±0.14 












1.0 


0.78 ± 0.43 








4>NR 






0.53 ±2.4 








-log£ 


86.0 


133.6 


83.6 







not improve the fit quality that much. The worst fit to the data is provided by a model with 
just a non-resonant amplitude. From this analysis, we estimate that the significance of the 
Z b (10610) -> BB* signal exceeds the 8a level. 

As the nominal model for the T( 10860) — > B*B*tt decay, we use the following parame- 
terization: 

S B *B*Tv{m) = (A Zb (w650) + A NR )E B * B * n (m). (4) 

We also fit the data without a non-resonant component and with a non-resonant amplitude 
alone. Results of the fits are shown in Fig. lU^b); numerical values are given in Table HI 

The best description of the B*B*n data is achieved in a model with only the Zb( 10650) 
amplitude included. The addition of a non-resonant amplitude does not provide any sig- 
nificant improvement of the fit quality. The fit with a non-resonant amplitude alone gives 
a much worse likelihood value. From this analysis, we determine the significance of the 
Z fe (10650) ->• B*B* signal to be 6.8a. 

In all fits discussed above, the masses and widths of the Zf, states were fixed at the 
values obtained from the analysis of the T(nS)7T + 7T~ and hb{mP)TT + 7T~ — final states: 
M[Z ft (10610)] = 10607.2 ±2.0 MeV/c 2 , r[Z 6 (10610)] = 18.4 ±2.4 MeV and M[Z b (10650)] = 
10652.2 ± 1.5 MeV/c 2 , r[Z 6 (10650)] = 11.5 ± 2.2 MeV. If allowed to float, the fit returns 
10597 ± 9 MeV for Z 6 (10610) mass in the fit to BB*ix events and 10649 ± 12 MeV for 
Zfc(10650) mass in the fit to B*B*n events. Large errors here reflect a strong negative 
correlation between resonance mass and its amplitude. 

ANALYSIS OF T(10860) -> T(n5)7r+7r- 

In addition to new results on the analysis of the three-body T(10860) — > B^B^n 
decays described in previous sections, we extend our analysis of the three-body T(10860) — > 
T(nS)7T + ir~ and T(10860) — > hb(mP)ir^ ~ir~ decays reported earlier in Ref. |jj to measure 
not only the parameters of the newly observed Zb states but also the fractions of individual 
components contributing to the three-body signals. 
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FIG. 5: Distribution of recoil mass associated with the ir + ir~ combination for T(5S) —> T(nS)ir + Tr~ 
candidate events in the (a) T(15); (b) T(2S); (c) T(35) mass region. Vertical lines define the 
corresponding signal region. 



To select T(10860) — > Y(nS') 7r+7r ~ ( n — 1, 2, 3) candidate events, we require the pres- 
ence of a pair of muon candidates with an invariant mass in the range of 8.0 GeV/c 2 < 
M(n + n~) < 11.0 GeV/c 2 and two pion candidates of opposite charge in the event. All 
tracks are required to originate from the vicinity of the interaction point. We also require 
that none of the four tracks be consistent with being an electron. More details on the 
analysis flow can be found in Ref. [lj and references therein. 

Candidate T(10860) — > Y(nS') 7r+7r ~ events are identified via the measured invariant mass 
of the combination and the recoil mass, M r (7r + 7r~), associated with the 7t + tc^ system 

calculated as M r (n + -K~) = ^/(E cm — E* +7c _) 2 — p* +„._, where E^+^- and p^+ n - are the 
energy and momentum of the 7r + 7r~ system measured in the cm. frame. Events originat- 
ing from T(10860) decays are selected with the requirement of \M r (-K + -n~) — M(n + p~)\ < 
0.2 GeV/c 2 . The M r (Tt + ix~) distributions shown in Fig. [5] are fit to the sum of a Crystal 
Ball function for the T(nS) signal and a linear function for the combinatorial background 
component. Results of the fits are shown in Fig. 

For the subsequent analysis, we select events around a respective T(nS) mass peak as 
shown in Fig. After all the selections are applied, we come up with 1819, 2219 and 588 
events for the T(lS')7r + 7r _ , T(2S')7r + 7r _ and T(3S)7t + tt~ final state, respectively. The frac- 
tions of signal events in the selected samples are determined from the fit to the corresponding 
M r {it + it~) spectrum. 

The amplitude analysis of three-body T(10860) —> T(nS)ii + 7Y^ {n = 1, 2, 3) decays is per- 
formed by means of an unbinned maximum likelihood fit. The distribution of background 
events is determined using events in the T(nS) mass sidebands. The variation of recon- 
struction efficiency over the phase space is determined using MC simulated signal events 
generated with a uniform distribution. 

We use the following parameterization of the T(10860) —> T(nS)n + n^ three-body decay 
amplitude: 

M{s u s 2 ) = Ai(si, s 2 ) + A 2 ( Sl , s 2 ) + A fo + A h + A NR , 

where sj = m 2 (Y(nS)ir + ), s 2 = m 2 (Y(nS)ii~). The amplitudes Ai and A 2 are S"-wave Breit- 
Wigner functions to account for the observed ^(10610) and -2^(10650) peaks, respectively. 
To account for the possibility for the T(10860) to decay to both Z + n~ and Z~n + channels, 
the amplitudes Ai and A 2 are symmetrized with respect to tt + and 7r~ interchange. Taking 
into account isospin symmetry, the resulting amplitude is written as 

A k = a k e i5k (BW{ Sl ,m k , T k ) + BW{s 2 , m k , T k )), 



11 




FIG. 6: Comparison of fit results (open histogram) with experimental data (points with error bars) 
for T(15)-7r + 7r _ events (left column), T(25)7r + 7r _ events (middle column) and T(3S)it + it~ events 
(right column) in the signal region. The hatched histograms show the background component. 



where the masses rrik and widths Tk (k = 1,2) are free parameters of the fit. Due to the 
very limited phase space available in T(10860) ->• T(nS)n + 7i- decays, T(10860) Z+tT 
and T( 10860) — > Z^n + amplitudes overlap significantly. We also include amplitudes Af 
and Af 2 to account for possible contributions from / (980) scalar and / 2 (1270) tensor states. 
We use a Breit-Wigner function to parameterize the /2(1270) and a Flatte function for the 
/o(980). The mass and width of the / 2 (1270) state are fixed at their world average values (8[; 
the mass and coupling constants of the / (980) state are fixed at the values denned from 



the analysis of B + ->■ K+ir+ir-; M(/o(980V) = 950 MeV/c 2 , g W7T = 0.23, g KK = 0.73 [10 . 

Following the suggestion given in Refs.jlll. fl2l |. the non- resonant amplitude A NR is pa- 
rameterized as 

A NR = a ■ e itf " + b ■ e^ r . S3) 

where S3 = m 2 (7r + 7r~), a nr , 6 nr , and 8% are free parameters of the fit (S3 is not an 
independent variable and can be expressed via s\ and s 2 but we prefer to keep it here for 
simplicity) . 

The logarithmic likelihood function L is then constructed as 

C = -2 M/dgStsi, s 2 ) + (1 - f sis )B( Sl , 82 )), 

where S(si, s 2 ) is formed from \M(si, s 2 )| 2 convolved with the detector resolution and / S i g is 
the fraction of signal events in the data sample. Results of fits to T(10860) — > T(nS)n + n^ 
signal events are shown in Fig. |6l where one- dimensional projections of the data and fits are 
presented. 

Results on the Zf, parameters are reported in Ref. [lj. Here, we report fractions of inter- 
mediate channels contributing to each three-body final state. The results are summarized in 
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TABLE II: Summary of results on fractions of individual quasi-two-body channels contributiong 
to Y(10860) — > T(nS)ir + iT~ three-body decays. 



Final state T(15)vr+^- T(25)tt+^- T(3S)tt + tt- 

S(^(1O61O)7T±)X0(^(1O61O) T(nS>±), % 2.54$™$}* 19-6±i;f±J : | 26.81™ ± 1.5 

B(ZT(10650)tt±)xB(^(10650) -> T(nS>±), % l-04±gjf ±g;g 5.77lJ;*£°;g ll.Q+£§±0.7 

B(T(nS)/ 2 (1270))xB(/ 2 (1270) ^ vr+vr-), % 15.6 ±1.4 ±2.1 2.81^;^°;| 

Total S-wave, % 89.2 ± 3.0 ± 2.4 105.6 ± 4.1 ± 2.6 45.6 ± 5.3 ± 0.8 



Table IIVt where the central values are determined from fits with the nominal model. Statis- 
tical uncertainties are determined from fits to multiple toy MC samples generated according 
to the nominal model. This allows us to account for correlations between various channels. 
In general, we find that all the three-body T(10860) — > T(nS)7T + ir~ decays are dominated 
by the .S-wave channels with some statistically significant .D-wave contribution. 

The dominant systematic uncertainty in the fractions of individual channels contributing 
to three-body decays comes from the model uncertainty. We estimate this uncertainty by 
fitting the signal with various modifications of the nominal model. For example, we vary 
the parameterization of the non-resonant amplitude or replace the T(riS , )/' 2 (1270) amplitude 
with a D-wave non- resonant component. 



ANALYSIS OF T(10860) -»■ h b (mP)ir+ir- 

In the analysis of the T(10860) — > hb(mP)n + 7T~ decays, we perform an inclusive re- 
construction of signal events utilizing the recoil mass, M r (7r + 7r~), associated with a 7r + 7r~ 
pair. The selection requirements are identical to those described in Ref. |2|. The continuum 
e + e _ — » qq (q = u, d, s, c) background is suppressed by a requirement on the ratio of the 
second to zeroth Fox- Wolfram moments P 2 < 0.3 [lij. We select 7r =1= candidates that orig- 
inate from the vicinity of the interaction point and are positively identified as pions based 
on the CDC (dE/dx), TOF and ACC information. We reject tracks that are identified as 
electrons. 

Because of the extremely high combinatorial background, a Dalitz analysis of the 
T(10860) — > hb(mP)7T + ir~ decay is challenging. Instead, we study the one- dimensional 
projection by fitting the M r (ir + TT') spectra in bins of the /^(lP)^ mass. We define the 
^(lP)^ mass as the recoil mass, M r (7r T ), associated with a single charged pion. We sym- 
metrize the distributions by combining the M r (7r + 7r~) spectra corresponding to M r {ir + ) and 
M r (7r~) bins and restrict the analysis to the M r (?r) > 10.40 GeV/c 2 {M r (ir) > 10.57 GeV/c 2 ) 
region for the hb{lP)ix + n~ (hi ) (2P)7r + 7T~) final state to avoid double counting. Each 
M r .(7r + 7r _ ) spectrum is fit to extract the /if, (IP) and the hb(2P) signal yields. The fit- 
ting function is a sum of a Crystal Ball function for the h\,(mP) signal and a Chebyshev 
polynomial for the combinatorial background. We also account for the T(2S) signal and a 
reflection from T(3S f ) — > T(1S)tt + tt~ decay. Details of this analysis can be found in Ref. (2)]. 
The hb(mP) yields as a function of the M r (ir T ) are shown in Fig. [7J where a clear two-peak 
structure is apparent for both /i6(lP)7r + 7r~ and /if,(2P)7r + 7r~ final states. 

We perform a \ 2 fit of the M r {jt) distributions to a coherent sum of two P-wave Breit- 
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FIG. 7: The yield (dots with error bars) of (a) hj,(lP) and (b) hb(2P), as a function of M r (7r). 
The solid histogram shows the results of the fit. 



Wigner amplitudes and a non-resonant contribution: 

\BW x {s, Mi, Ti) + ae^BWi(s, M 2) T 2 ) + 6e # | 2 ^, (5) 

v s 

where \fs = M r (7r); the variables M&, T^. (A; = 1, 2), a, 0, 6 and ^ are free parameters; ^ is a 
phase-space factor, where p (q) is the momentum of the pion originating from the T( 10860) 
(Zt>) decay measured in the rest frame of the corresponding mother particle. The P-wave 

Breit- Wigner amplitude is expressed as BWi(s, M, T) = where F is the P-wave 

Blatt-Weisskopf form factor F = yj ^j^p - (lH| . go is the daughter momentum calculated 

using the pole mass of its mother, R = 1.6 GeV -1 . The function in Eq. (j3J) is convolved with 
the detector resolution function (er = 5.2 MeV/c 2 ), integrated over the histogram bin and 
corrected for the reconstruction efficiency. The result of the fit is shown in Fig. [7J fractions 
of individual amplitudes are given in Table IIHI 

The non-resonant contribution in the hb(lP)ix + ii~ final state is found to be consistent 
with zero, b = 0.181q56 with the central value for its fit fraction of only 3.2% and an upper 
limit of 22% at 95% C.L. This is in accordance with the expectation that the non-resonant 
amplitude is suppressed due to the heavy quark spin flip. In the hb{lP)Ti + Ti~ final state, 
we are not sensitive to the non-resonant amplitude due to the limited phase space and thus 
fixed the non- resonant component at zero. 

To estimate the systematic uncertainty we vary the order of the Chebyshev polynomial 
in the fits to the M r (7r + 7r~) spectra. To study the effect of finite M r (7r) binning, we shift 
the binning by half of the bin size. To study the model uncertainty in the fits to the 
M r (7r) distributions, we remove (add) the non-resonant contribution in the hb(lP) (/t&(2P)) 
case. We also vary the R parameter of the Blatt-Weisskopf form factor in the range of 
— 5 GeV -1 (default value is 1.6 GeV -1 ) and find the associated systematic effect to be 
negligible. In Ref. 0], we find that the resolution in data could be larger than in MC by 
typically 10%; we increase the width of the resolution function by 10% to account for possible 
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TABLE III: Fit fractions of various components in three-body T(10860) —> hb(rnP)Tr + n signals. 





h b (lP)ir+ir- 


h b (2P)ir+ir- 


non-resonant 
Z 6 (10610)7r ± 
Z 6 (10650)7r ± 


3.2% (< 22% at 90% C.L.) 
(42.3i? 2 5 7 ^; 8 7 )% 
(60.2±^;f ±| : |)% 


(35.2™ 
(64.8^ 1? 5 7 5 )% 



difference between data and MC simulation. The maximum change of parameters for each 
source is used as an estimate of its associated systematic error. All systematic uncertainty 
contributions have been added in quadrature to obtain the total systematic uncertainty. 



RESULTS 

To calculate branching fractions for the observed three-body T( 10860) — > BB*ir and 
T(10860) — > B*B*7r signals, one needs to account for the non-uniform distribution of signal 
events over the phase space. The corrected efficiencies are found to be 12.25 ± 0.06% and 
11.0 ± 0.1%, for the T(10860) ->■ BB*n and T(10860) -> B*B*n modes, respectively. The 
three-body branching fractions are then calculated as 



i3(T(10860) B^B^n) 



N B M B (*) 7 



L ■ a(e+e- -> T(10860)) ■ B(B -> /) • b bWbWv • a 



where L = 121.4 fb _1 is the total integrated luminosity and a(e + e~ — > T(10860)) = 0.340 ± 
0.016 nb (9| is the cross section of e + e~ annihilation into the T(10860) state at a cm. energy 
of 10865 MeV. Using world average results for the secondary branching fractions |8J|, the 
combined fraction of B meson decays to all reconstructed final states including secondary 
branching fractions is found to be BiB —>•/) = (143 ± 15) x 10~ 5 (neutral and charged B 
combined). Finally, one also needs to correct for the oscillation of neutral B mesons. After 
time integration, the fraction of oscillated neutral B mesons is equal to / Q sc = 0.19. The 
correction factor a is calculated as: 

B(B+ -> /+) ■ e BWflWir + B(B° -> f) ■ e flWflWir (l - / osc ) 
a = - 



B(B+ -> /+) • e BMBMir + B(B° /°) • eflWBW, ' 

where B(B + — > f + ) and B(B° — > f°) are the total fractions of B decays to charged and 
neutral final states (including secondary fractions), respectively. From signal MC, one gets 
a = 0.8978. This results in i3(T(10860) BB*n) = (28.3±2.9±4.6) x 10" 3 and T(10860) -»■ 
B*B*tc = (14.1 ± 1.9 ± 2.4) x 10^ 3 . For the T(10860) ->■ BBn decay, we calculate a 
90% confidence level upper limit of T( 10860) — > BBtt < 4.0 x 10~ 3 (including systematic 
uncertainty). 

The dominant sources of systematic uncertainties for the three-body branching frac- 
tions are the uncertainties in the secondary branching fractions, the uncertainty in the 
reconstruction efficiency and in the signal yield extraction and the uncertainty in the 
a(e + e~ — > T(10860)) cross section. The overall systematic uncertainties for the three-body 
branching fraction are estimated to be 17.5%, 16.3% and 16.9% for the BB7r, BB*7r and 
B*B*tt final states, respectively. 
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TABLE IV: Results on three-body T(10860) — > T(nS)ir + Tr branching fractions. 

Final state T(lS)ir + ir- T(25)vr+7r- T(35)vr+^- 

Signal Yield 2090 ± 73 2476 ± 97 628 ± 41 

Efficiency, % 45.9 39.0 24.4 

B(T(nS) -)• [8], % 2.48 ±0.05 1.93 ±0.17 2.18 ±0.21 

B(T(10860) -> /), 10~ 3 4.45 ±0.16 ±0.35 7.97 ± 0.31 ± 0.96 2.88 ± 0.19 ± 0.36 

From Ref. [13], 10" 3 5.3 ± 0.3 ± 0.5 7.8 ±0.6 ±1.1 4.8+^ ± 0.7 



The branching fractions of the three-body T(10860) — > T(nS)7i + 7i decays are calculated 
with the following formula: 

£(T(10860) -> T(nS)7i + 7i-" - r(nS)n+n- 



L ■ a(e+e- -> T(10860)) • B(T(nS) -> • e T (nS)»+ w - ' 

The reconstruction efficiencies (including trigger efficiency) £-r(nS)7r+7r- are determined from 
the signal MC, T(nS) — > fractions [8|; the final results are given in Table IIVI In 

determination of the reconstruction efficiencies, we use signal MC events generated ac- 
cording to the results of the Dalitz fit with the nominal model. The main systematic 
uncertainties in three-body fractions come from the a(e + e~ — > T(10860)) cross section 
(4.7% for all channels), the T(nS) — > jjL + n^ branching fractions (2.0%, 8.8% and 9.6% 
for n = 1,2,3, respectively), the T(nS) signal yield (4.5%, 5.3% and 4.9% for n = 1,2,3, 
respectively), and the MC tracking efficiency of 4% for all channels. The overall systematic 
uncertainty is 7.9%/12.0%/12.4% for T(lS/2S/3S)n + n- , respectively. The results for the 
T(10860) — > T(nS)7r + 7T~ fractions are to be compared with previous measurements by Belle 
with a data sample of 21 fb _1 [13| (see Table HV|) . We find the two sets of measurements to 
be consistent within uncertainties. 

Using results of the fit to the M r (7r) spectra with the nominal model (see Tabled]) and 
the results of the analysis of the T(10860) -»■ T(nS), n = 1,2,3 and T(10860) -»■ h b {mP), 



TABLE V: List of branching fractions for the Z^(10610) and -2^(10650) decays. 



Channel 


Fraction, % 






Z 6 (10610) 


Z 6 (10650) 


T(15)vr+ 


0.32 ±0.09 


0.24 ±0.07 


T(25)vr+ 


4.38 ± 1.21 


2.40 ±0.63 


T(35)vr+ 


2.15 ±0.56 


1.64 ±0.40 


h b (lP)TT+ 


2.81 ± 1.10 


7.43 ± 2.70 


h b (2P)7T+ 


4.34 ±2.07 


14.8 ±6.22 


B+B*° + B°B*+ 


86.0 ±3.6 








73.4 ± 7.0 
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m = 1, 2 decays, one can measure the ratio of the branching fractions: 



6(^6(10610) BB*) 



6.2 ± 0.7 ± 1.3 



+0.0 
-1.8 



£ n £(Z b (10610) ->• T(n5)7r) + ^ m Z 6 (10610) -)• h b (mP) 



and 



£(Z fe (10650) ->■ 5*5*) 



2.8 ±0.4 ±0.6 



+0.0 
-0.4- 



E n B(Z 6 (10650) T(nS)Tr) + £ m Z 6 (10650) ^(mP) 



We also find it useful to calculate the relative fractions for Z b decays assuming that thy 
are saturated by the already observed T(nS) (n = 1,2,3), hb(mP) (m = 1,2), and B*B^ 
channels. The results are summarized in Table IV! We do not include the ^(10650) — > BB* 
channel in the table as this decay mode has marginal significance. However, if the central 
value is used, its fraction would be 25.4 ± 10.2%. All other fractions would be reduced by a 
factor of 1.33. 
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